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ABSTRACT 


A survey of the inner Galaxy region of Galactic longitude / G [+15°, +50°] 
and latitude b G [—4°, +4°] is performed using one-third of the High Altitude 
Water Gherenkov (HAWG) Observatory operated during its construction phase. 
To address the ambiguities arising from unresolved sources in the data, we use a 
maximum likelihood technique to identify point source candidates. Ten sources 
and candidate sources are identihed in this analysis. Eight of these are associated 
with known TeV sources but not all have differential fluxes compatible with 
previous measurements. Three sources are detected with signihcances > 5 u after 
accounting for statistical trials, and are associated with known TeV sources. 

Subject headings: astroparticle physics - gamma rays; diffuse background - gamma 
rays: general 
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Introduction 


We present a TeV gamma-ray survey of the inner Galaxy region between -f 15° and -|-50° 
in Galactic longitude (/) and ±4° in Galactic latitude (6) using data collected with a partially 
completed conhguration of the High Altitude Water Gherenkov (HAWG) Observatory. This 
region of the sky contains the strongest Galactic detections in the held of view of HAWG 
other than the Grab Nebula, and is known to have a large number of gamma-ray sources 
such as pulsars, pulsar wind nebulae (PWNe), supernova remnants (SNRs), compact object 
binaries, diffuse emission from the Galactic plane, and sources without known astrophysical 
associations. The most commonly identihed TeV Galactic sources are SNRs and PWNe 


(Wakely & Horan 2008). SNRs are thought to be acceleration sites of Galactic cosmic rays. 


since they can accelerate particles via diffusive shock acceleration and provide sufficient 


power to explain the cosmic-ray energy losses from the Galaxy (e.g. Drury et ah (2001)). 
In PWNe, electrons effectively gain energy at the termination shock where the pulsar wind 
is terminated by the surrounding gas, emitting TeV gamma rays via inverse Gompton 


scattering (Aharonian et ah 2008a). Many of the gamma-ray sources in the Galactic plane 


are unidentihed (HID) due to the fact that the measurements can envelope multiple sources 
identihed at other wavelengths. Morphological and spectral studies are crucial for making 
associations with observations at other wavelengths and for distinguishing leptonic and 
hadronic gamma-ray production processes which will aid in source identihcation. 

The HAWG survey presented here overlaps in both energy range and sky coverage with 
the survey and observations performed by the H.E.S.S. imaging atmospheric Gherenkov 


telescopes (lAGTs) (for example, Aharonian et al. (2006b, 2008b); [Hoppe (2008); Djannati- 


Atai et al. (|2008); Ghaves et al. (2008[)). Over 15 sources within this region were discovered 


by H.E.S.S. Similar surveys within this region have been performed by Milagro (Abdo et 


al. 2007) and ARGO (Bartoli et al. 2013b), along with targeted observations by the lAGTs 
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VERITAS and MAGIC (for example, Albert et al. (2006); Acciari et al.[(2009); Aleksic et 


al. (2012); Alin et al. (2014a); Aleksic et al. (2014)). At lower energies, the Fermi Large 


Area Telescope (FermALAT) has pnblished its third sonrce catalog in the lOOMeV-300 GeV 


energy range (3FGL) based on the hrst fonr years of science operation (Acero et al. 2015). 
The region snrveyed with the partial HAWC array inclndes the locations of 73 sonrces 
from this catalog, 47 of which are withont known astronomical associations. The hrst 
catalog with sonrces > 10 GeV (IFHL) has been pnblished based on the hrst three years of 
FermALAT data (Ackermann et al. [2013 ). The region snrveyed in this pnblication contains 
twelve sonrces from this catalog, fonr of which are not in the 3FGL catalog. 

Sonrce identihcation is a major challenge when analyzing the emission from this region, 
with point-like and extended emissions overlapping each other. For this reason, a method is 
applied to the data to simnltaneonsly ht the positions and diherential hnx normalizations of 
mnltiple sonrces assnming a simple power law with a spectral index of 2.3. The maximnm 


likelihood method used here, similar to the one used by the FermALAT (The Fermi-LAT 


Collab. 2012), is described in more detail in Section 3, following a description of the HAWC 


detector in Section 2. In Section 4, we present a list of 10 sources and candidate sources 
with locations and differential flux normalizations found with the likelihood method, 
followed by a discussion of their possible associations with previously reported objects. 


2. The HAWC Gamma-Ray Observatory 


The HAWC Gamma-Ray Observatory is located at Sierra Negra, Mexico (18°59’41” 
N 97°18’30.6” W) at 4,100 m a.s.h, and is sensitive to gamma rays and cosmic rays in the 


energy range from 100 GeV to 100 TeV (Abeysekara et al. 2013a). The observatory consists 


of an array of 300 water Cherenkov detectors (WCDs) as shown in Fig.[^ covering an area 
of 22,000m^. Each WCD consists of a tank 7.3 m in diameter, 4.5m in depth, and hlled 
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with 230,000 L of purified water. Four upward-facing photomultiplier tubes (PMTs) are 
attached to the bottom of each WCD: one high-quantum efficiency 10-inch Hamamatsu 
R7081-MOD PMT at the center and three 8-inch Hamamatsu R5912 PMTs spaced 120° 


apart 1.8 m from the center (Abeysekara et ah 2013b). The array has an instantaneous field 
of view of 2 sr and a duty cycle > 95% is expected once construction is completed. 


The PMTs detect Cherenkov light inside the WCDs produced by relativistic secondary 
particles from an extensive air shower. The PMT signals propagate through 600 feet of 
RG59 cable to the data acquisition system in a central building, where they are processed 
by custom-made front end boards. The PMT pulses are shaped and discriminated using a 
low and a high voltage threshold. The time stamp of a pulse that crosses at least the lower 
of the two thresholds is recorded by CAEN VX1190A time-to-digital-converter modules 
with a precision of 0.1 ns. The pulse size and ultimately the charge can be inferred from the 
time-over-threshold (ToT). 


A laser calibration system measures the ToT-charge conversion and the response time 


of each PMT to different pulse sizes (Abeysekara et ah 2013b). A simple multiplicity 
trigger requiring at least 16 to 21 PMTs above threshold within a 150 ns time window is 
used to identify air shower events. Once an air shower triggers the detector, the charge and 
timing calibration are applied to each PMT signal. The Nishimura-Kamata-Greisen (NKG) 


function (Greisen 1956 Kamata & Nishimura 1958) describes the lateral distribution of 
charged particles produced by an electromagnetic air shower. The distribution of signal 
amplitude in the detector is fitted with the NKG function in order to find the position of 
the shower core, where the extrapolated trajectory of the primary particle intersects the 
altitude of the detector. The direction of the primary particle is subsequently reconstructed 
by taking the estimated position of the shower core and fitting the relative times of the 
PMT signals produced when the air shower crosses the array. 
















x[m] 


Fig. 1.— The layout of the full HAWC array, with filled circles marking the PMTs operational 
at the beginning of the dataset under investigation and small open circles indicating the 
PMTs that were eventually integrated into the array during this period. WCDs are shown 
as large circles. WCDs not commissioned during the period covered by this dataset are 
drawn without PMTs. 
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3. Data and Analysis 

The HAWC array was completed in March 2015, but science operations began in 
August 2013. The analysis described in this paper is performed on data taken with the 
partially constructed HAWC array (called Pass 1 data/conhguration hereafter) between 
August 2, 2013 and July 9, 2014. During this period, data taking was occasionally 
interrupted for detector construction work and maintenance, and the duty cycle is ~ 84%. 
The dataset contains 275 ± 1 source transits of the inner Galaxy region, depending on 
the right ascension. The detector grew from 362 PMTs in 108 WCDs to 491 PMTs in 
134 WCDs during this time as shown in Fig.[^ The data analyzed are divided into three 
epochs, distinguished by the number of active PMTs. 

The sensitivity of HAWC is a function of source declination. The best sensitivity is 
achieved for sources that transit through the detector zenith. For the study presented here, 
a zenith angle cut of 45° is applied, corresponding to declinations between —26° and +64°. 

The energy, angular resolution, and background rejection of the observed gamma rays 
are correlated with the shower size measured in the array. Therefore, the data are divided 
into 10 bins according to the fraction / of PMTs triggered by an air shower event passing 
standard selection cuts (e.g. > 1 PE within 450 ns time window) out of the total number of 
active PMTs in the array. Cuts are applied to three parameters in each / bin to separate 
the cosmic-ray background from the gamma-ray signal: 

1. The ratio of the chi-square of the £t to the core location and the number of PMTs in 
the event; 

2. A topological cut based on the “compactness” of the charge distribution in the array; 

3. The reduced chi-square of the £t to the shower direction. 
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The cuts are explained in detail in Abeysekara et al. (2015a). 


The primary background of this analysis is air showers induced by hadrons. The shape 
of the lateral distribution of gamma-induced air showers differs from that of hadron-induced 
showers. The NKG function used for fitting the shower core describes the lateral 
distributions of gamma-ray showers. Consequently, the chi-square values resulting of the 
core £t are on average smaller for gamma-ray showers than for cosmic-ray showers with the 
same number of triggered PMTs, and this parameter is used for background discrimination. 
Furthermore, hadronic showers produce pions with large transverse momenta that decay 
to gamma rays and muons. These subshowers are likely to produce large signals in PMTs 
far from the shower core. Fig.[^ shows a typical gamma-like event and a hadron-like event 
observed with the Pass 1 configuration. The ratio of the number of triggered PMTs to the 
number of photoelectrons (PE) in the PMT that detects the strongest signal outside of a 
radius of 40 m from the shower core is found to result in good gamma/hadron separation 
performance. A third selection applied to the reduced chi-square distribution of the shower 
angle fit removes poorly reconstructed air showers. All cuts are optimized by maximizing 
the sensitivity of the Pass 1 data to emission from the Crab Nebula, the brightest steady 
source at TeV. In the lowest / bin, about 79% of gamma ray events and 38% of cosmic ray 
events pass the cuts, whereas for the highest / bin these numbers drop to about 13% and 
0.03% respectively. 


3.1. Sky Maps and Maximum Likelihood Method 


The arrival directions of air showers that pass the cuts are binned in a signal map in 


equatorial coordinates using HEALPIX (Corski et al. 2005), a pixelization scheme that 


splits the unit sphere into twelve equal-area tessellations. Each of the twelve tessellations 
are subdivided into an N^ide x N^ide grid, giving 12 N^^de pixels. In this analysis, N^ide is 
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Fig. 2.— A gamma-like event (left) and a hadron-like event (right) observed by the Pass 1 
conhguration. Each hlled circle is a PMT signal, with the color indicating the hit time and 
the size indicating the charge. The reconstructed shower core is marked with a red star and 
a dashed circle indicates the 40 m radius around the reconstructed core. The highlighted red 
circle represents the location of the PMT with the maximum observed PE outside of the 
40 m radius from the core. 
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chosen to be 512, which provides an average pixel size of 0.11°. The bin size is mnch smaller 
than the typical point response of the Pass 1 detector, discussed in Section 3.2. 


The background at each location in the sky is estimated from the data using the direct 


integration method described in Atkins et ah (2003). It is computed by the convolution of 
the local event hour angle and declination distribution with the all-sky event rate recorded 
during a predehned integration step. An integration step of two hours is used in this 
analysis, corresponding to 30° in right ascension, to emphasize structures smaller than this 
angular scale. 


To convert the gamma ray counts to flux, a simulation-based model of the detector 
response to sources of gamma rays is used. The detector response hie describes the energy 
distribution and the point spread function (PSF) as a function of source declination and 
fraction hit / bins. It is generated with the air shower simulation program CORSIKA (|Heck 


et al.j 

1998 

) and the detector simulation package GEANT4 (Agostinelli et ah 

2003) 


the data are subdivided into three epochs rehecting the size of the Pass 1 array, three sets 
of detector simulations with 365, 426, and 480 PMTs were generated. 

The hrst step of the maximum likelihood ht is to build a source model characterized 
by the position and spectrum of the source. Right ascension a and declination 6 are used 
to describe the source position. The source spectrum in the present analysis is assumed to 
obey a simple power law 


dN _ 
dE ~ ° 



( 1 ) 


where Jq is the differential hux normalization, Eq is the pivot energy, which is chosen at 
where the differential flux normalization is least dependent on the spectral index, and P is 
the spectral index. The study presented here assumes a hxed index of 2.3 due to the limited 
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sensitivity in this data set. The index of 2.3 is representative of measured values for known 
Galactic objects. 

The region of interest (ROI) used for a likelihood £t needs to be larger than the angular 
resolution of the detector to include most photons from a given source. However, it is not 
always possible to hud an ROI with photons from one single source due to a high potential 
of source confusion in the Pass 1 Galactic plane data. Therefore, the source model may 
need to contain more than one source and in this case the expected event count becomes 

^ij = bij + ^ 7p’fc; (2) 

k 

where bij is the background events in the jth pixel of the Rh / bin, and ■jijk is the expected 
number of gamma rays from the kth source in the jth pixel of Rh / bin. The event count 
is convolved with the detector response. As the observed event count in each pixel is 
distributed according to a Poisson distribution, the probability of observing N number of 
events given an expected count A from the source model is 

A^e-^ 

WA) = ^^. (3) 

The likelihood given a parameter set 9 = {a, 6, Iq) in the source model is the product 
of the likelihood of each pixel in an ROI and in each / bin: 


f bins ROI 

£(?iiv)= n w 

i j 

where Nij and Xij are the observed and expected event count in the jth pixel of the ith / 
bin, respectively. The logarithm of the likelihood is used for ease of computation: 
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fbins ROI 


\nCie\N) = - A,,), 


( 5 ) 


The term Nij\ is discarded from Eq.[^ since it is independent of the parameters in the source 
model. The log likelihood is maximized with respect to the parameter set 6 in the source 
model using the MINUIT package (|Brun fc Rademakers 1997). 


A likelihood ratio test is performed to decide how many sources are needed to properly 
model an ROI. To decide if the one-source model is preferred over the background-only 
model, the log likelihood of the background-only model ln£o is computed hrst. Then the 
log likelihood In Ci of the one-source model is computed. The test statistic (T S) defined by 


TS = -2{\n Co-In C^) (6) 

is used to compare the goodness of £t between the two models. The T S'-value is converted 
to a p-value, which is the probability of the data being consistent with the background-only 
hypothesis. The same likelihood ratio test can be used to compare between two models 
with N and N-l-1 sources. 

In this iterative process, an additional source characterized by three free parameters 
(right ascension, declination, and differential flux normalization) is added to the model 
at each step, with the positions and amplitudes of the existing sources free to change in 
response to the new source. We repeat this procedure as long as ATS > 15 after adding 
a new source, corresponding to a p-value of 1%. The ATS threshold was chosen a priori. 
After the iteration in each ROI, a simultaneous £t of the differential source fluxes in all 
regions is performed while the source positions are fixed. This is done to take into account 
photons from a source with centroid position just outside of an ROI but still contributing 
to the ROI due to the PSF of the detector. Finally, the differential flux and TS value of 
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each source are obtained by fitting a single source while treating other sources as part of 
the background with hxed positions and differential fluxes. 

Fig-i shows the pre-trials signihcance map of the inner Galaxy region surveyed in this 
paper. The map is made by moving a putative point source through each pixel, performing 
a maximum likelihood £t of differential flux normalization with the spectral index fixed at 
2.3, and converting the TS value to a signihcance according to Wilks’ theorem (a = y/TS 
with 1 degree of freedom (DoF) for this map). There are multiple > 5cr hotspots in the 
surveyed region as shown in Fig.[^ The area is divided into hve ROIs as dehned in Table[2 
In each of these ROIs, the procedure described above is applied to simultaneously account 
for the hux contributions of neighboring sources. 


3.2. Systematic Uncertainty and Simulation Studies 

The energy range in which the present analysis is most sensitive is estimated using 
detector and air shower simulations. Since more inclined air showers travel through more 
atmosphere, the median gamma-ray energy of the Pass 1 data is a function of declination. 
The declination dependence is shown in Fig.for the spectral index assumption F = 2.3. 
The median energy increases from 7TeV for a declination of -1-19° to ~ 30TeV for 
declinations of —26° and -|-64° assuming a spectral index of 2.3. 

The Crab Nebula is used to measure the point spread function in the data since its 
extent is much smaller than the PSF of the HAWC detector. We describe the PSF of the 
Pass 1 detector as the linear combination of two Gaussian functions. Fitting the emission 
from the Crab Nebula region in each / bin, the 68% containment region is found to vary 
from 2.5° to 0.6° depending on the event size that correlates with the fractional number of 
PMTs in the event (see Fig.j^on 68% PSF as a function of / bin). The PSF measured in 
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Fig. 3.— Top: Significance map; Middle: Model with 11 seed sources (crosses) and no 
uniform surface brightness £t. Black circles indicate the 5 regions of interest; Bottom: 


Residual signihcance map. The open squares mark TeV sources in Wakely & Horan (2008). 
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Table 1. Definition of the 5 regions of interest 


Region 

RA 

Dec 

Radius 

1 

286.4° 

6.7° 

3.0° 

2 

284.2° 

2.5° 

3.0° 

3 

280.7° 

O 

1 

O 

o 

4 

279.6° 

-9.2° 

o 

O 

5 

275.9° 

-13.1° 

3.0° 



Fig. 4.— Median gamma-ray energy of the Pass 1 data as a function of declination along 
with the 10% and 90% quantiles, derived from simulations and assuming a spectral index of 


-2.3. 
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the data from the Crab Nebula is used for all other sources in this likelihood analysis. 


The pointing of the detector is verihed using data from the Crab Nebula and the 
blazar Markarian (Mrk) 421. Both source regions are relatively isolated in the sky and do 
not suffer overlap with additional sources. A likelihood £t with right ascension, declination, 
and differential flux normalization is performed using a power law spectral assumption 
with the hxed spectral index of 2.6 for the Crab Nebula and 3.0 for Mrk 421. The results 
are summarized in Table|^ The Crab Nebula spectral index assumption is chosen based 
on previous lACT measurements, and for Mrk 421 a soft spectral index of 3.0 is chosen 
due to known spectral cutoff that is not modelled here. The positions of both sources in 


the Pass 1 data are consistent with measurements by lACTs (for example Aharonian at 


al. (2004); Albert et ah (2007)). Changing the spectral index used in the £t of the Crab 
Nebula between 2.0 and 3.0 shifts the best £t position by < 0.07°. The signihcance maps in 
the vicinity of the Crab Nebula and Mrk 421 are shown in Fig.|^ These maps are made by 
moving a putative point source through each pixel and performing a maximum likelihood £t 
of the differential flux normalization with spectral index hxed at 2.6 and 3.0, respectively. 
Then the TS value in each pixel is converted to signihcance according to Wilks’ theorem. 


The Crab Nebula is the most significant source in the Pass 1 data. The spectrum of 


Table 2. Pointing in J2000 


Source 

SIMBAD Database 

RA (°) Dec (°) 

RA"^ (°) 

Pass 1 

Dec*^ (°) 

TS 

signihcance 

Crab 

83.63 

22.01 

83.53 ±0.06 

22.06 ±0.06 

491.4 

22.2 a 

Mrk 421 

166.11 

38.21 

166.22 ±0.18 

38.14 ±0.18 

69.0 

7.8 a 
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Fig. 5.— PSF measured in data from the Crab Nebula (68% containment) as a function 
of the fraction / of PMTs used in event reconstruction. The systematic uncertainty on the 
measured PSF (gray) is < 20%. A linear combination of two Gaussians is needed to correctly 
characterize the PSF. 
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Fig. 6.— Significance maps of the Crab Nebnla (left) and Mrk 421 (right). The white 


sqnares indicate the source positions from Wakely & Horan (2008) and the white crosses are 


measured positions in the Pass 1 data along with 1 a errors. 


this source has been well measured by lACTs. The Pass 1 data cover a similar gamma-ray 
energy range as the lACT data. The pivot energy for the Crab Nebula analysis is selected 
to be 4 TeV in order to minimize the dependence of the differential flux normalization on 
the spectral index in the £t. It differs from the pivot energy of ITeV used by Whipple, 
HEGRA, H.E.S.S., and VERITAS and 0.3 TeV used by MAGIC, so the lACT differential 
flux normalizations at 4 TeV are computed from the respective flux normalizations and 
spectral indices ( Hillas et al.||199^ [Aharonian at al. 2004 Aharonian et al. 2006a Alin et 
al. 2014b Albert et al. 2008[ ). Tab le[^ shows the Pass 1 differential flux normalization from 
the Crab Nebula with a spectral index assumption of 2.3, which is used for the analysis of 
the Inner Galaxy region, as well as the differential flux normalization with a spectral index 
assumption of 2.6, which is close to the index measured by lAGTs. The differential flux of 
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the Crab Nebula as measured in the Pass 1 data is within 15% of the lACT measurements. 
The flux of Mrk 421 is beyond the scope of this paper due to known variability of the blazar 
during the period covered by this dataset and will be presented in detail in an upcoming 
publication. 

The determination of the differential flux normalization in Pass 1 depends on the 
assumption of the spectral index, but the variance can be minimized by selecting the best 
pivot energy. For example, the derived Crab differential flux normalization increases by 5% 
if an index of 2.3 is chosen instead of 2.6 for a pivot energy of 4TeV, while a pivot energy of 
9 TeV results in a 30% increase in flux normalization. For the sources in the region of the 
Galactic plane analyzed here a spectral index of 2.3 is assumed and the best pivot energy 
varies from 4 TeV to 9 TeV depending on declination. In order to minimize the variance 
in the flux normalization determination due to different spectral indices, different pivot 
energies are used for sources as a function of declination. According to simulations, this 
leads to a systematic change of —10% to +20% in the htted differential flux for a source 
spectral index ranging between 2.1 and 2.8. 

There are three other major contributions to the systematic error of the flux 
normalization that have been studied using data and simulations: 

1. Detector conhguration variability. The number of active PMTs changed continuously 
in this dataset but only three conhgurations are modeled with the detector 
simulations. The event passing rate is found to vary < 20% among the three simulated 
conhgurations, which is equivalent to the resulting uncertainty on the hux estimate. 

2. Angular resolution. The measured PSF on the Crab Nebula is used to compute the 
expected number of gamma rays in each pixel in a given model. The error in the 
measured PSF width from the Crab Nebula is < 20%, which corresponds to 15-20% 
uncertainty in the hux estimate. 
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Table 3. Differential Flux Normalization Comparison of the Crab Nebula Assuming a 

Simple Power Law 


Instrument 

Flux Normalization at 4TeV 

( 10-13 TeV“^cm“^s“^) 

Reported Flux Normalization 

at 1 TeV^ (10“^^ TeV“^cm“^s“^) 

Spectral Index®' 

HAWC Pass 1 

8.57 ±0.45 

_ 

2.30'" 

HAWC Pass 1 

8.25 ±0.40 

- 

2.60'" 

Whipple 

10.1 

3.2 ±0.17 ±0.6 

2.49 ±0.06 ±0.04 

HEGRA 

7.49 

2.83 ±0.04 ±0.6 

2.62 ±0.02 ±0.05 

H.E.S.S. 

9.00 

3.45 ±0.05 ±0.7 

2.63 ±0.01 ±0.09 

MAGIC 

9.25 

57 ± 2 ± 6" 

2.48 ±0.03 ±0.2 

VERITAS 

8.83 

3.48 ± 0.14 ± 1.08 

2.65 ±0.04 ±0.3 


'^The first quoted uncertainty is statistical uncertainty and the second is systematic uncertainty. 
’^Assumed spectral index. 


At 300 GeV. 
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3. Charge scale. The core htter and gamma/hadron cuts rely on the charge measurements 
by PMTs, which are based on the ToT-charge calibration. Up to 20% change in 
the flux estimate is observed due to the uncertainty in the charge scale estimated in 
studies of single muons with HAWC. 


Each source of systematic uncertainties contributes ~ 20%. In addition, a minor 
contribution of ~ 8% from the uncertainty of atmospheric modeling is taken into account 
(Abeysekara et ah 2015b| ). We add systematic uncertainties in quadrature for a total 
systematic uncertainty in the flux normalization of ~ 40%. 


4. Results 


Table|^ lists the epoch J2000 positions, differential flux normalizations, TS, and 
the post-trials significances of the detections and candidates from this analysis. Eleven 
seed sources are initially identified with the ATS > 15 criterion and are used in the 
source model. To calculate the number of trials within the inner Galaxy search region 
of 280 square degrees, a Monte Carlo study of fluctuations based on the measured 
background was performed to obtain the p-value distribution of the maximum significance 
on background-only maps within the inner Galaxy region. A fit to the corresponding tail 
distribution has been performed, and the exponent derived from this fit is equal to the 
number of trials. The resulting number of trials is 424 ± 3. After accounting for trials 
, ten source detections and candidates remain with > 3 a. Table[^ lists the 
possible counterparts of each source and differential flux normalization comparison with 
known TeV sources. The discussion is separated into TeV source detections and source 
candidates below using a criterion of 5 a significance after trials. 


(Biller 1996 










Table 4. Galactic Source Detections and Candidates 
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Table 5. Possible TeV Gamma-Ray Source Counterparts 
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After the identification of 11 seed sonrces, a model and a nniform snrface brightness 
for the entire region was htted simnltaneonsly to the data. While for the sonrces a spectral 
index of 2.3 is assnmed, the spectral index assnmption for the nniform snrface brightness is 
2.5. The ATS of adding the nniform snrface brightness as another free parameter to the 
source model is 33, i.e. 5.7a that the uniform surface brightness component is preferred. 
The htted surface brightness at 5TeV is (1.6 ± 0.4) x 10“^^ TeV“^cm“^s“^sr“^, which is 
compatible with the average diffuse hux of (1.0 ±0.2) x 10“^^ TeV“^cm“^s“^sr“^ reported by 


H.E.S.S. extrapolated to 5TeV within the same region (Abramowski et al. 2014). However, 
the uniform surface brightness measured in this dataset is not simply diffuse emission but 
also a combination of unidentihed sources, source extensions, and photon contaminations 
from sources due to uncertainties in the PSF. As evident in the residual map of Fig|^ there 
are several 3 a regions around known TeV sources that are not detected in this dataset and 
are contributing to this uniform surface brightness ht. The uniform surface brightness ht 
is not included in TablelH The contribution to the source differential hux normalization is 
< 30% of the smallest reported hux. 


4.1. Source Detections 

The source IHWC J1857±023 is detected at 6.2 a post trials and is ~ 0.4° away from 
both HFSS J1857±026 and HESS J1858±020. These two TeV sources were discovered by 
the H.E.S.S. collaboration during their Galactic plane survey and are ~ 0.7° apart. The hux 
of HESS J1857±026 is approximately an order of magnitude higher than HESS J1858±020 


(Aharonian et al. 2008b), and the diherential hux normalization from IHWC J1857±023 


is compatible with the combined hux of HESS J1857±026 and HESS J1858±020. Both of 
the HESS sources were detected as extended, with HESS J1857±026 as the larger of the 


two. Aleksic et al. (2014) reported energy dependent morphology for HESS J1857±026 with 











two distinct components, MAGIC J1857.2+0263 and MAGIC J1857.6+0297. These sources 
cannot be resolved with these data from the partial HAWC array. The spectrum reported 
by MAGIC is for the entire region and is compatible with the differential flux normalization 
derived from this dataset. 


IHWC J1838-060 is detected at 6.1 a post trials and is located in the middle of the 
known TeV sources HESS J1837-069 and HESS J1841-055. This detection overlaps with 
the extension of HESS J1841-055, and the differential flux normalization is compatible with 


that reported by H.E.S.S.(Aharonian et ah 2008b). Bartoli et ah (2013a) reported a 5.3a 
detection by ARGO near the position of HESS J1841-055, with a 0.4° source extent and is 
closest to this detection. The ARGO source extends towards HESS J1837-069 and includes 
several Ferm+LAT sources. The flux reported by ARGO, when converted to differential 
flux at 7TeV for comparison, is ~ 4x the differential flux normalization derived from this 
dataset. 


IHWC J1825-133 has a post-trials signihcance of 5.4 a. It is ~ 0.5° to the south of the 
HESS J1825-137 centroid position, which is an extended PWN with spectral softening as a 


function of distance from the pulsar towards a southeast direction (Aharonian et ah 2006c). 
The simple power-law flux derived from this dataset is lower than the flux extrapolated 


from the simple power-law assumption measured by H.E.S.S. However, Aharonian et ah 


(2006c) reported the spectrum is unlikely to be a simple power law and presented several 
alternative hts. The derived flux normalization from this dataset is most compatible with 
the energy dependent photon index power law £t by H.E.S.S. There is also a nearby PSR 


0.4° away, PSR J1826-1256, seen by Ferm+LAT (Acero et ah 2015) and associated with 


the Eel nebula (Ferrand & Sah-Harb 


2012 ). 
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4.2. Source Candidates 


IHWC J1907+062c is 4.6 a post trials with a best-fit position that is compatible with 


previously reported positions of MGRO J1908-I-06 (see Abdo et ah (2007); Aharonian et 


al. (2009); Bartoli et ah (2012); Alin et al. (2014a) for example). The differential flux 
normalization is consistent with the flux measured by H.E.S.S. and VERITAS and in 
agreement with Milagro given the statistical uncertainties of both instruments. Alin et ah 


(2014a) reported strong excess near the pulsar PSR J1907-I-0602 but also extends toward 


SNR G 40.5-0.5. The Pass 1 dataset is not able to resolve the spatial morphology of this 
source. 

IHWG J1904-|-080c has a post-trials significance of 3.9 cr. There is currently no 
previously reported TeV detection near this location. The nearest gamma-ray source is 


3FGL J1904.9-1-0818 at 0.3° away (Acero et al. 2015). However, this is a weak detection 


from the Fermz-LAT 3FGL catalog, at < 5 a, with no known association. 

IHWG J1844-031c has a post-trials signihcance of 4.7 a and is spatially coincident 


with HESS J1843-033, which is classihed as an unidentified source (Hoppe 2008). However, 


the morphology of this detection appears to extend towards HESS J1846-029, a pulsar 


wind nebula (Djannati-Atai et al. 2008). Bartoli et al. (2013b) reported a 4.2cr excess. 


ARGO J1841-0332, associated with HESS J1843-033 despite being 0.7° away, due to the 
large systematic pointing error at high zenith angle. 

IHWG J1849-017c is detected at 3.7 a post trials and is positionally coincident with 
the extended source HESS J1848-018, which is possibly associated with the star forming 


region W43 (Ghaves et al. 2008). The differential flux normalization at 6 TeV from this 


dataset is ~ 3.5x the flux reported by H.E.S.S. Ghaves et al. (2008) reported an index of 


2.8 and the spectral index assumption of 2.3 in this analysis would result in a different flux 
normalization by 20%. More importantly, diffuse emission from this star forming region 
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that contains a molecular cloud could contribute more to the differential ffux normalization 
derived from the Pass 1 dataset than that measured by lACTs, which have a smaller 
angular integration region. 

IHWC J1842-046c has a post-trials significance of 3.4 a and has no clear gamma-ray 


association. There is an X-ray source, SNR G27.4-1-0.0 (Kes 73), located ~ 0.4° (Ferrand 


& Safi-Harb 2012) away from this candidate. 


IHWC J1836-090c is detected at 3.9 a post trials. It is spatially coincident with 
HESS J1834-087 and the SNR W41 ( Aharonian et al.]|2006b Albert et al.||2006 Abramowski 
et al.|2015). The differential flux normalization from the Pass 1 dataset at 8TeV is ~ 6x 


higher than the flux reported in Abramowski et ah (2015). The source is reported by 


H.E.S.S. as having a central point-like component and an extended component. A similarly 
extended component is also seen by FermALAT. The region contains a candidate pulsar at 


the center of the SNR W41, and two scenarios are supported by Abramowski et ah (2015): 


PWN or SNR interaction with a nearby molecular cloud. The cloud density traced by 
appears wider than the H.E.S.S. detection and may contribute to the increased flux seen in 
the Pass 1 data with HAWC. 

IHWC J1836-074c has a post-trials significance of 3.2 a, with the nearest TeV PWN, 


HESS J1837-069 (Aharonian et ah 2006b), ~ 0.5° away with a compatible differential flux. 


There is also a GeV source 3FGL J1837.6-0717 (Acero et ah 2015) that is ~ 0.3° away with 
no association. 


5. Conclusion 


A survey of the inner Galaxy has been presented in the region of Galactic longitude 
I G [15°,50°] and latitude b G [—4°,-|-4°] using 283 days of data from the partial HAWC 
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Gamma-Ray Observatory from August 2013 to July 2014. Three sources have been detected 
at > 5 a with an additional seven candidate sources detected at > 3 a after accounting for 
trials. While associations with previously published lACT detections are not always within 
the experimental uncertainties, about half of them have differential flux normalizations that 
are compatible with the previous detections. 

A likelihood method similar to the Fermi-LAT source-hnding algorithm has been 
applied to data from an extended air shower array for the hrst time to properly address 
challenges arising from source identihcation. The point sources presented here have 
differential fluxes > 20% of the Crab Nebula flux at several TeV. Future analyses will build 
on this method to most effectively exploit the increased sensitivity and improved pointing 
of data collected with the full HAWC array, and to explore other regions of the sky visible 
to the observatory. 
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